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Abstract 

In the present paper we discuss the hmits and the correct utihzation of the 
standard program for the inversion of Ehashberg equations and the determina- 
tion of the electron-phonon spectral function and the coulomb pseudopotential 
from tunneling measurements in high-Tc superconductors. In order to com- 
pare the calculated density of states with the experimental one, we introduce 
the results of the inversion procedure, applied to our recent tunneling data in 
Bi2Sr2CaCu208+x single-crystal break junctions with Tc = 93 K, in a direct 
program for the solution of Eliashberg equations. Most of the observed dif- 
ferences between theoretical and experimental curves at energy greater than 
gap can be explained by a smooth energy dependence of the normal density of 
states that we introduced in the direct solution of the Eliashberg equations. 
Finally we show that the effects of the energy-dependent normal density of 
states can be simulated by an efficient electron-phonon spectral function but, 

also, by a negative, nonphysical, coulomb pseudopotential. 
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I. INTRODUCTION 



At the present time it can be regarded as a well-established fact that in the oxide su- 
perconductors, as in conventional ones, and at temperatures below the critical temperature 
Tc, the charge carriers are bound in pairs in an energy band whose width is of the order of 
the gap, but the mechanism responsible for the attraction between the charges has still not 
been identified. 

A very important method for obtaining information concerning this aspect is constituted 
by tunnel spectroscopy measurements |T|-0, that, through the Eliashberg equations (EE) 
inversion, allow to determine the electron-phonon spectral function, o? (cu) F (cu), of the ma- 
terial. Seven years after the publication of the first, tentative, o?[ijj)F{uS) curves obtained 
by inversion of the gap equations from point-contact tunneling data in Bi2Sr2CaCu208+3; 
(BSCCO) [0,^, a certain skepticism still persists on the possibility of obtaining a reliable 
and reproducible o?[lS)F{uS) from tunneling data in BSCCO. The main reasons for this skep- 
ticism are related to: i) the determination of the low-temperature normal-state conductance 
that is a fundamental ingredient for the inversion of the gap equations, but it is not a measur- 
able quantity in this high-Tc superconductor and ii) the considerable deviations of the dl /dV 
curves from the ideal BCS behaviour, observed in this material up to short time ago. More 
recently, tunneling experiments using the break-junction approach have shown that single 
crystals of very high structural and crystallographic quality and a very fine technique for 
the production of stable and reproducible junctions can lead to BCS-like lifetime-broadened 
results In this case the deviations of the BSCCO break-junction dl/dV characteris- 

tics from the ideal Superconductor-lnsulator-Superconductor (SIS) behaviour are relatively 
small. After an unfolding procedure for extracting the density of states Ns from the SIS 
tunneling conductance, it is possible to observe that most of the deviations of Ns from the 
BCS density of states can be explained in all the energy range only by the presence of a 
certain amount of lifetime broadening 0. Similar results have been recently obtained with 
the point-contact technique by using highly doped GaAs as counterelectrode [Q]. 
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In this work we develop a method for the correct utihzation of the best experimental tun- 
neling data in the inversion of the EE, by introducing different theoretical energy-dependent 
normal densities of states (NDOS) determined by suitable criteria we will discuss in detail in 
section III. In this way we have calculated a'^{uj)F{ijj), the coulomb pseudopotential ji* and 
other observable physical quantities that depend on the first two, for every normal density 
of states. Then we have put a'^{uj)F{uj) and ii* in a direct program for the solution of the 
EE and, for each different NDOS, we calculated the superconductive quasiparticle density 
of states (DOS) for comparison with the experimental data. Some possible explanations 
for the discrepancies between theoretical and experimental DOS are given. By using the 
direct program and including the effect of an energy-dependent symmetric NDOS of proper 
shape, we find that the agreement with the experimental data improves very much at ener- 
gies greater than gap. Finally, from the theoretical DOS calculated by the direct program 
with an energy-dependent NDOS, and by using again the standard inversion program we 
determine the efficient values of a^{uj)F{ijj) and that simulate the effects of the NDOS 
depending on energy. From now on we set ^ = c = /c^ = 1. 



II. INVERSION OF THE ELIASHBERG EQUATIONS 

The experimental data that we use in order to extract the a^{(jj)F{oj) and to compare 
with the theoretical DOS have been obtained by employing the break- junction technique on 
BSCCO single crystals of very high crystallographic and chemical quahty, as indicated by 
x-ray diffractograms, resistivity versus temperature and susceptibility measurements. Before 
breaking the crystals, their R{T) characteristics were measured and the critical temperature 
was Tc = 93 K with Arc(10 - 90%) = 2 K. Immediately afterwards, at 4.2 K, a finely con- 
trolled bending force was applied to the thin sample holder which the crystal was mounted 
on. It produced the break of the very thin crystal (about 1 x 1 x 0.02 mm^) perpendicularly 
to the ah planes and the creation of the junction. By varying the bending force it was 
possible to reproducibly and accurately control the resistance of the junction and also the 
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I{V) curves were very stable and reproducible. All the details concerning the experimental 
data and the reproducible determination of the density of states and the Eliashberg function 
in these crystals can be found in a forthcoming paper [^]. Among several I{V) tunneling 
measurements of very high quality we selected one of the best and applied an unfolding 
procedure for extracting the quasiparticle density of states Ns from the SIS tunneling con- 
ductance. Figure 1 shows that the quality of the experimental DOS (open circles) is largely 
improved with respect to the measurements of few years ago and, in particular, sharp peaks 
at the energy gap A and an almost constant Ns at eV ^ A are now present. Nevertheless, 
it is possible to observe that deviations of Ns from the "pure" BCS density of states Nbcs 
are still present: it is a problem for the inversion of EE p|-|l4| because very large oscillations 
in Nrid = Ns/ {Nn ■ Nbcs) — 1 (where N^ is the normal density of states) produce singularity 
and no convergence at the increase of the number of iteractions when the standard inversion 
program is used [0,0]. Usually this difficulty in the calculation of the reduced density of 



states Nrid is overcome by using the lifetime-broadened expression [3-p|JH] instead of Nbcs- 



Ni{uj) = real ^{lu — iV)/ \J {uJ — iVy — A(0)^| . The values of the energy gap A(0) and of 
the lifetime broadening parameter F that fit the experimental density of states of Fig. 1 are 
A(0) = 23 meV, F = 2.9 meV and, consequently, 2A{0)/T^ = 5.74. In Fig. 1 the solid line 
represents the lifetime broadening fit of the experimental data. 

Even if widely used as a first approximation for the a'^{u!)F{uj) determination, the pro- 
cedure of using Ni{uj) instead of Nbcs appears somehow incoherent. In fact, some parts of 
the inversion program contain the ideal reference curve, Nbcs, but, for calculating N^id, a 
different curve A'; is used. For a fully coherent utilization of the inversion program it is nec- 
essary to compare the experimental data to the ideal BCS behaviour. We also believe that 
it is fundamental to normalize the experimental data to an energy- dependent NDOS and 
then calculate the Nrid by using the Nbcs as reference. In high-T^ superconductors (HTS) 
like BSCCO, the normal conductance is unknown at the temperature of measurement, for 
example Tea, = 4.2 K, since the upper critical magnetic field necessary to suppress the su- 
perconducting phase is very large and unrealizable, while, on the other hand, AT = Tc — T^x 



is too large to consider NDOS(Tc)^ ND0S(Te2:). For this reason, we will use four simple 
arbitrary expressions for the NDOS with three free parameters each, to be determined by 
suitable physical conditions. Comparing the results of this approach to the experimental 
data of Fig. 1 we will select the NDOS that, possibly, will be the most plausible one for our 
break-junction tunneling experiment. 



III. ENERGY DEPENDENCE OF THE NORMAL DENSITY OF STATES 



The HTS have a very complicated crystallographic structure. Ho et al. ||T5[ have given 
general motivations to the fact that six or more atoms of the transition elements for unitary 
cell may produce definite peaks in the NDOS. There are also many experimental results 
indicating this hypothesis as correct [^,|TB|,|T^. In particular the results of Mandrus et al. 
in BSCCO [|l^ and our recent measurements in Bi2Sr2Cu06+2; single crystals suggest an 
energy dependence of the NDOS which varies with the direction of the tunneling process 
{ab plane or c axis). Since the exact dependence is unknown, we used four different, but 
symmetric, analytical forms for the energy dependence of the NDOS: exponential, linear, 
parabolic and lorentzian 
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where N{0), a and x are free parameters. A^(0) is the density of states at the Fermi level and 
X controls the behaviour of the DOS at large energies. In the best tunneling measurements 
along the ab plane, it is reasonable to consider x > 0, i.e. the tunneling conductance is 



constant or slightly decreasing at large energies, both for experimental reasons P,^,|8|jr8| and 
because the total isotope effect is less than 1/2 in HTS and this implies d'^Nniui) / doj"^ U=o< 
||T9| . Of course in all the Eqs. 3.1-3.4, A^„(ti;) = A^(0) when uj tends to zero. 

For determining the values of the parameters we employ the following conditions: 

a) \\m^^+^Ns{uj) = \im^^+^ N^iuj) = Nexi+oo) 

where Nex{+oo) is the experimental density of states at energy u ^ A(0); 

b) /o+°° N,x{u)du; = /o+~ Nn{u)du 

which expresses the law of conservation for the number of states; 

c) Jl^){NUu^)/[N4uj) ■ Nbcs{u^)] - l}du - 

since the oscillations of the experimental "normalized" data around the ideal BCS curve 
in the phonon energy range should be symmetric |ll|,pO|. 

The condition b is ideal: in reality the experimental curve shows a small zero-bias con- 
ductance and we have applied this modified condition: 

b') /o"'" N,x{u)du - Ai, = J^- N^{uj)duj 

where is the highest energy of the experimental data and Ai^ = c ■ Um with < c < 
Nex{0) is a small contribution eventually produced by leakage currents. 

In the same way, for the practical estimation of the parameters, the condition c becomes: 
c') min \jl^{N,x{uj)/[N4uj) ■ Nbcs{uj)] - l}duj 

where > A(0) is the first point of the experimental data used in the standard 
procedure for the inversion of the Eliashberg equations. 

In the case of our experimental data shown in Fig. 1, a program for the numerical 
calculation of the parameters according to the conditions a, b' and c' yields : 





iV(0)(fi-i) 


X 


a^meV) 


Ai,{meV ■ fi-i) 


exp 


6.099 ■ 10-2 


0.108 


52 


0.403 


lin 


6.635 • 10-2 


0.206 


70 


0.288 


par 


6.230 ■ 10-2 


0.132 


69 


0.345 


lor 


6.200 ■ 10-2 


0.127 


35 


0.368 
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The different, energy-dependent NDOS curves corresponding to the parameters of the 
table are shown in Fig. 2. If we know the I{V) curve for T > Tc, in order to check the 
quahty of the fit for Nn{uj), we can use the well-known expression of the tunnel current of 
a symmetric junction: 

HV) = e-'G^N ■ N^{uj)N^{uJ + eV) [f{uj) - f{u + eV)] du 

where Gnn is a constant and f{uj) is the Fermi function. Since Gnn is unknown, we 
can use the quantity I{V)/[dI{V)/dV] for a comparison of the theoretical results with the 
experimental data. 

Now we can correctly calculate the Nrid = Ng^i^^) /[Nn{uj) ■ Nbcs{^^)] — '^ and, by the stan- 
dard inversion program [|l^, we can determine the spectral function a"^ {u) F (u) , the coulomb 
pseudopotential /i*, the electron-phonon interaction constant A = jQ°°[a'^{Lj)F{uj)/uj]du!, T^, 
the area A of a'^{uj)F{u), the coefficient b of the quadratic part of a'^{uj)F{u!) for — 
and Uc, that is the cut-off frequency. The results of this inversion procedure, applied to 
the experimental data of Fig. 1, are summarized in the following table for the different 
energy-dependent NDOS curves: 





A 




A{meV) 


Tc{K) 


b{meV 2) 


ojcijneV) 


exp 


3.34 


0.93 ■ 10-2 


47.21 


93.6 


3.9-10-3 


225 


lin 


2.86 


0.26 ■ 10-1 


47.68 


99.7 


1.6 ■ 10-3 


225 


par 


3.03 


0.31 ■ 10-1 


48.40 


98.8 


2.7-10-3 


225 


lor 


3.24 


0.31 ■ 10-1 


47.49 


94.6 


3.0- 10-3 


225 



In Fig. 3 the spectral functions a^{uj)F{uj) for the different NDOS models are compared 



to the generalized phonon density of states G{uj) of BSCCO pT |. 

By solving the standard Eliashberg equations in direct way we can now check if the 
different a'^{uj)F{uj) and /i* reproduce the experimental density of states Nex{oj). The well- 
known standard EE are: 



Piuj)- 



(3.5) 



7 



Jo 



■oo 



[1 - Z{uj)]-uj 



(3.6) 



where 




uj'+uj + Q + i- 6 uj'— uj + Q — i-6 



and P{uj) = Real {/S.{ijj) / Joj"^ — /S?{ijj)\is the pair density of states, while Ns{uj) = 



Figures 4 (a), 4 (b), 5 (a) and 5 (b) show a comparison of the densities of states cal- 
culated for Nn{uj) = cost by solving Eqs. 3.5 and 3.6 (dash-dot lines) with the normalized 
experimental ones (open circles) and with the densities of states calculated by the modified 
EE that take into account the four different models for Nn{uj) as it will be discussed in the 
following (solid lines). Discrepancies between the theoretical DOS for Nn{uj) = cost and the 
normalized experimental curve are remarkable for all the different forms of Nn{uj). They 
can be due to the inversion procedure as well as to the form of the Eliashberg equations. 

For what concerns the inversion program, it is necessary to remember that: 

1) The true normal conductance can be more cumbersome than the theoretical one and, 
as a consequence, Nnd may be not correct. 

2) The standard inversion program was made for ideal BCS-like curves, not for gapless, 
broadened ones. 

3) The dispersion relations, used in the standard program for calculating Real{A{u!)) 
and Im{A{uj)) are incorrect if iVex(O) ^ 0. 

4) In the standard inversion program lim^^Q a"^ {u) F (u) = lim^^o^ ' (parabolic be- 
haviour) whereas in superconductors with bidimensional character like BSCCO we can ex- 
pect a linear behaviour: lim^^o (y'^{uj)F{u!) = lim^^o d ■ OJ. 

5) The energy dependence of the normal density of states is neglected in the standard 
inversion program, i.e. Nn{uj) = N{0). 

Concerning the Eliashberg equations it should be born in mind that: 
1) The Migdal's theorem is no more valid in HTS since the Debye phonon frequency is 
of the order of the Fermi energy |2^ . 



Real yuj/Juj'^ — A^(c<j))is the quasiparticle one. 
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2) The Eliashberg equations in the simple form expressed by Eqs. 3.5 and 3.6 are vahd 
only for homogeneous and isotropic three-dimensional superconductors whereas BSCCO is 
bidimensional and very anisotropic. 

3) The order parameter is supposed to be in s-wave symmetry. 

4) The coulomb pseudopotential is considered constant but it could be energy dependent 
2|. 

5) In the simple, standard form of the Eliashberg equations Nn{uj) = N{0). 
In spite of the previous assertions, we assume two fundamental hypotheses: 

a) The standard inversion program provides results approximately correct for what con- 
cerns a'^{uj)F{u!) and /i*. 

b) The differences between experimental and theoretical curves at u; > A can be ex- 
plained by putting Nn^uj) ^ A^(0). 

The new Eliashberg equations that take into account the energy dependence of the normal 



density of states at T = are p^- p^] : 



l\{oj)Z{uj) = / 77-1 (u; ) ■ Kj^ioj^ijj^ — li* du , 

Jo L -I 

[1 — Z(uj)] ■ CO = / n2{uj ) ■ K^{uj,uj )duo , 
Jo 

Jo 

The normal density of states depending on energy p(u;) = Nn{u!)/N{0) enters through 
three functions defined as: 



n2{uj) = { --) Im( f p{uj') ■ j^T^^du' 

71/ \J-oo L>{UJ, iO 



naiuj) = (-) Im( f p{uj') ■ ^^^-p^dJ 

\7T/ \J~oo D[UJ,UJ ) ^ 



where D{lu,uj') = Z\uj) ■ [uj^ - A\uj)] - [lu' + xi^)f ■ 

Energies are measured from the chemical potential and the real part of x(^) represents a 
shift of the chemical potential due to the electron-phonon interaction. The energy- dependent 
electronic density of states p^u) modulates the pair density of states A^(0) ■ ni{u) and the 
quasiparticle one A^(0) ■ [^2(1^) — n-^^uj)]. In this article, only for simplicity, we have supposed 
that all the normal densities of states are symmetric: Nn{uj) = Nn{—uj) and this produces 
X(^) = 0. 



IV. RESULTS AND DISCUSSION 

The theoretical Ns{uj)/N{0) curves for Nn{uj) 7^ cost (solid lines) are compared to the 
experimental normalized curves (open circles) in Figs. 4 (a, b) and 5 (a, b) for exponential, 
linear, parabolic and lorentzian cases, respectively. The same figures also show the normal 
densities of states A^„(ci;) (dot lines) and Nbcs{^)/ + (dash lines) for the four above 
mentioned models. We can see that the agreement between the theoretical and the experi- 
mental normalized DOS is particularly good for the exponential and the linear cases (Figs. 
4 (a) and 4 (b), respectively). 

The gap value and the ratio 2A(0)/Tc are modified by the energy dependence of the 
NDOS as it is shown in the following table: 





A{0){meV) 


2A(0)/T, 


exp 


21.0 


5.24 


lin 


20.0 


4.99 


par 


21.2 


5.27 


lor 


21.3 


5.31 



As it has been predicted in a previous work a negative curvature of the energy- 
dependent NDOS produces a reduction of the gap A(0). 

The best case seems to be the exponential one where A = 3.34, /i* = 0.93 ■ 10^^, 
Tc = 93.6 K. We remind that the experimental value of Tc is 93 K. In our opinion, 
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these results are remarkable because fi* is small but positive, Tc is coincident with the 
experimental one and A is not too large: the amorphous low-T^ superconductor Pbg.sBio.s 
has A = 3.00 and 2A(0)/Tc = 5.19 [110]. Moreover if we calculate the electron-phonon 



coupling constant by using the following self-consistent formula p3| , |29| for T = K: 
A = 2 /q^°° doj'Nniyj') /q^°° dVL "^^J*^^)^'' obtain A = 2.69. It is correct to remember that Tc 
is calculated by the EE with Nn{oj) = N{0), but, since Nn{uj) has not a definite peak near 
the Fermi level, the value of Tc is not significantly affected by the energy dependence of the 
NDOS [^. In Fig. 6 the experimental data (open circles) are compared to the theoretical 
electron density of states for the exponential case (solid line) in the full positive energy 
range. The agreement of the two curves for uj > 25 meV is very good. 

Similar results have been obtained in 6 additional break-junction tunneling curves mea- 
sured on the same BSCCO single crystals with Tc = 93 K. The dl/dV curves suitable for the 
analysis previously described were selected among all the measured data by using some cri- 
teria that can be summarized as follows: i) Presence of sharp and symmetric gap structures; 
ii) Presence of symmetric and clear phonon-like structures; iii) Presence of a background 
conductance smoothly decreasing with energy at > 25 meV. In all the cases the use of 
the exponential energy dependence of Eq. 3.1 for the normal density of states gave the best 
agreement between the theoretical DOS and the experimental one. The electron-phonon 
spectral functions a'^{u!)F{u!) obtained by the inversion of the Eliashberg equations in these 
6 additional cases show large similarities to the a'^{uj)F{uj) of Fig. 3 (solid line) and the 
calculated average values of A and Tc present a small variance of the order of ±11% and 
±2.7%, respectively g. 

An interesting question is whether it is possible to reproduce the tunneling characteristics 
for Nn{uj) 7^ cost within the framework of the usual Eliashberg theory. For trying to give 
an answer we considered the Nrid{uj) calculated by the direct program for Nn{uj) = N^^^{u) 
and, by using the standard inversion program, found an effective electron-phonon spectral 
function a"^ (u) F (u) cj f and a that simulate the energy dependence of the normal den- 
sity of states. The results are shown in Fig. 7 where the effective electron-phonon spectral 
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function is represented as a dash line while the original a'^{uj)F{uj) determined from ex- 
perimental data (after normalization with N^^^{uj)) is shown as a solid line and the open 
circles represent the neutron generalized phonon density of states. We obtained Ae// = 3.4, 
= -0.31 • 10~\ Tceff = 94 K, Aeff = 40.7 meV, befj = 9 • lO'^mel^-^ 

Finally, Figure 8 shows the electron-phonon coupling factor a'^{uj) as determined from the 
data shown in Fig. 7. We can observe that a'^{uj) determined from experimental data (solid 
line) has a large energy dependence which suggests the presence of a strong electron-phonon 
coupling with phonon modes between 12 and 27 meV and between 60 and 80 meV, while the 
coupling appears rather depressed for phonon modes between 29 and 56 meV. Similar results 
have been obtained in previous break-junction tunneling experiments in BSCCO and can be 
explained by the possible directionality mainly along the ab plane of the tunneling process in 
our break junctions Moreover, in Fig. 8 it can be seen that not to consider the energy 
dependence of the NDOS (i.e. to use the standard program for the inversion of the EE on the 
DOS obtained by the direct solution in presence of the NDOS energy dependence) produces 
an effective coupling factor a'^{uj)eff (dash line) that overestimates the soft phonon modes 
and gives rise to a little and negative nonphysical value of the coulomb pseudopotential. 

In conclusion, we have seen that in order to obtain plausible values for a'^{uj)F{uj), A, fi* 
and Tc it is necessary to use the standard program for the inversion of EE in a correct and 
coherent way. The effect of a nonconstant NDOS on the normalized tunneling conductance 
is significant and cannot be accounted for in a meaningful way within the framework of 
the standard Eliashberg theory without obtaining effective parameters a'^{uj)F{uj)eff and 
/igjj that alter the physical interpretation of the problem. Finally, a strong electron-phonon 
interaction, restricted to some particular phonon modes of the G{uj), and a smooth energy 
dependence of the NDOS permit to explain our recent break-junction tunneling data in 
BSCCO at > A in an excellent way. In the next future, it will be necessary to include 
in the EE the effect of anisotropy and of the bidimensional character of Bi-based high-Tc 
superconductors, as well as, to calculate the effect of the breakdown of Migdal's theorem 
on the critical temperature and on the superconducting gap by using the electron-phonon 
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spectral functions a^{u!)F{u!) reproducibly determined in this and in previous experiments. 
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FIGURES 

FIG. 1. Quasiparticle density of states determined from the tunneling conductance of a sin- 
gle-crystal Bi2Sr2CaCu208+a: break junction. 

FIG. 2. Comparison between the four different energy-dependent normal densities of states 
used in this paper and described in Eqs. 3.1 - 3.4. 

FIG. 3. Spectral functions a^{u)F{oj) determined by using the standard program of inversion 
of the Eliashberg equations for the different NDOS models of Fig. 2 compared to the generalized 
phonon density of states G{io) (open circles). 

FIG. 4. (a) Comparison of the calculated DOS for Nn{u) = cost (dash-dot line) with the 
experimental one determined from the data of Fig. 1 (open circles) and with the density of states 
calculated according to the NDOS exponential model of Eq. 3.1 (solid line). The BCS density of 
states (dash line) and the NDOS (dot line) are also shown for completness; (b) The same as in (a) 
but for the linear model of the NDOS. 

FIG. 5. (a) The same as in Figs. 4 (a) and 4 (b) but for the parabolic model of the NDOS; (b) 
The same as in (a) but for the lorentzian model of the NDOS. 

FIG. 6. Experimental DOS of BSCCO (open circles) compared to the theoretical electron 
density of states for the exponential NDOS model of Eq. 3.1 (solid line). 

FIG. 7. Comparison of the effective electron-phonon spectral function (dash line) to the original 

a'^{Lo)F{uj) determined from experimental data (solid line) and the neutron generalized phonon 
density of states (open circles). 

FIG. 8. Comparison of the original a'^{u) determined from experimental data in the frame- 
work of exponential NDOS model (solid line) and the effective electron-phonon coupling function 

a^(a;)e// (dash line). 
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